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We present a study of phase transition to macroscopic superfluidity for an ultracold bosonic gas 
confined in a combined trap formed by a one-dimensional optical lattice and a harmonic potential, 
focusing on the critical temperature of this system and the interference patterns of the Bose gas 
released from the combined trap. Based on a semiclassical energy spectrum, we develop an analytic 
approximation for the critical temperature T c , and compare the analytic results with that obtained 
by numerical computations. For finite temperatures below T c , we calculate the interference patterns 
for both the normal gas and the superfluid gas. The total interference pattern shows a feature 
of "peak-on-a-peak." As a comparison, we also present the experimentally observed interference 
patterns of 87 Rb atoms released from a one-dimensional optical lattice system in accord with our 
theoretical model. Our observations are consistent with the theoretical results. 
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I. INTRODUCTION 

Bosonic atoms confined in optical lattices have proved 
to be a unique laboratory for investigating quantum 
phase transitions from superfluids to Mott insulators [J- 
3;]. The momentum distribution of a lattice system can 
be mapped out directly by the interference pattern of 
the atomic cloud after a ballistic expansion over a time 
of flight (TOF). The emergence of macroscopic bosonic 
superfluid is usually identified by the appearance of in- 
terference peaks. However, recent theoretical works 
for homogeneous gases in a three-dimensional (3D) lat- 
tice showed that this criterion of macroscopic superflu- 
idity is not reliable since even a normal gas can have 
sharp interference peaks. The underlying physical pic- 
ture is that a lattice system at finite temperatures pos- 
sesses a "V-shaped" phase diagram [4-6] which includes 
a normal gas region between the Mott Insulator and the 
superfluid. The true signature of macroscopic superflu- 
idity is the (5-function momentum peaks with nearly unit 
visibility [4j. Below critical temperature, the coexistence 
of superfluid and normal gas in the homogeneous lattice 
system should give rise to an interference pattern hav- 
ing a feature of "peak on a peak" @ . The new criterion 
of macroscopic superfluidity makes it necessary to fur- 
ther investigate the phase transition of bosonic atoms in 
an optical lattice, particularly for the characteristics as- 
sociated with the critical temperature and interference 
patterns. Experimental investigations are also required 
for comparison with relevant theoretical models. 

There have been a few theoretical works consider- 
ing the translationally invariant (uniform) lattices. How- 
ever, in a realistic experiment, an optical lattice is always 



accompanied by harmonic confinement in all dimensions, 
arising from the focused Gaussian laser beams and/or an 
external magnetic trap. A bosonic gas is, therefore, never 
spatially uniform over the lattice range. Wild et al. (Toj 
have examined the critical temperature of the interacting 
bosons in a one-dimensional (ID) lattice with additional 



harmonic confinement. Ramakumar et al. [ll| have in- 
vestigated the condensate fraction and specific heat of 
non- interacting bosons in ID, two-dimensional (2D), and 
3D lattices in the presence of harmonic potentials. Based 
on a piecewise analytic density of states extended to ex- 
cited bands, Blakie et al. [12j developed an analytical ex- 
pression of the critical temperature for an ideal bosonic 
gas in the combined harmonic lattice potential, and com- 
pared the analytic result with their numerical compu- 
tations. However, these studies on combined traps did 
not mention interference patterns of the released bosonic 
gases. A more recent theoretical paper [l3[ has investi- 
gated the Bose-Einstein condensation (BEC) in a 3D in- 
homogeneous optical lattice system, and predicted that a 
bimodal structure in the momentum-space density profile 
is a universal indicator of BEC transition. 

The experiment of Spiclman et al. 



14| has examined 
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the superfluid to normal transition for a finite-sized 2D 
optical lattice system. Their measurements confirm that 
bimodal momentum distributions are associated with the 
superfluid phase. For such a system with a typical density 
of 1 atom per lattice site, the phase transition behavior 
can be interpreted in terms of the commonly used Bose- 
Hubbard model. 

Unlike the 2D and 3D cases, an inhomogeneous ID 
optical lattice system is usually much more heavily pop- 
ulated, with an atomic number up to several hundreds 
in a single lattice site. In the superfluid phase, the on- 
site interaction energy U varies from site to site because 
of its dependence on the local population in single lat- 
tice sites. This increases the complexity in searching for 
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an analytical description of the phase transition. In this 
paper we present a study of the critical temperature and 
interference patterns of an ultracold bosonic gas confined 
by a ID optical lattice and an additional magnetic po- 
tential. The interference patterns of the normal gas and 
the condensed atoms are treated separately. The super- 
position of the two parts gives rise to a feature of "peak 
on a peak." Different from a homogeneous lattice system, 
however, the normal gas can-not produce sharp interfer- 
ence peaks. Furthermore, the theoretical results are com- 
pared with our preliminary experiment for a ID lattice 
system of 87 Rb Bose-Einstein condensates. 

Our theoretical model relevant to the phase transition 
is for ideal bosonic gases. In fact, interatomic interaction 
may affect the shape of the interference pattern, espe- 
cially for the condensed part which has a higher atomic 
density. In order to obtain a better match with the exper- 
iment result, we take the interaction energy into consid- 
eration for the condensed atoms during the time of flight. 
The computed result shows that interference peaks aris- 
ing from the condensed atoms can be significantly broad- 
ened due to the interaction effect. 

This paper is organized as follows. In Sec. HH we be- 
gin with a semiclassical energy spectrum for a combined 
harmonic lattice trap. Under the tight-binding approxi- 
mation and in the low energy limit, we derive an analyt- 
ical expression of the critical temperature for the atoms 
condensed to a superfluid state. The accuracy of the an- 
alytical results are checked with respect to the numerical 
calculations. Section Mil gives a description on how the 
interference patterns are calculated for the normal gas, 
as well as the Bose-condensed gas. In Sec. IIV1 we briefly 
introduce the experiment, and present the observed in- 
terference patterns for a comparison with our theoretical 
results. Finally in Sec. [Vj we summarize the obtained 
results. 



II. CRITICAL TEMPERATURE 

We now consider an ideal Bose gas confined in a 3D 
harmonic potential with axial symmetry around the z 
direction. The axial and transverse trapping frequen- 
cies are lu z and uj x = io y = u>±, respectively. Moreover, 
we assume that the axial confinement is much weaker 
than the radial confinement (cj 2 <C wi), so that the Bose 
gas is made cigar-shaped. A ID optical lattice potential, 
Vo sin 2 (fcz), is applied along the z axis, where k — ir/d is 
the wave vector of the lattice light, d denotes the lattice 
period, and Vq denotes the potential depth of the lattice. 
Vq can be written in terms of the recoil energy E r , say, 
Vo = sE r , where E r = h 2 k 2 /2m, and m is the atomic 
mass. The harmonic potential, together with the optical 
lattice, forms a combined trap written as 



In practice, an optical lattice is usually produced by a 
retro-reflected Gaussian laser beam which also produces 
a transverse confining potential, that can be simply ab- 
sorbed into w_l if it is non-negligible. 

To obtain the eigenenergies of the combined trap sys- 
tem, one needs to derive the single-particle Hamiltonian 
of the system and then numerically diagonalize it fl2l ]. 
Despite its accuracy for ideal Bose gases, this numer- 
ical method can not provide an analytic expression of 
the energy levels. The energy spectrum corresponding 
to the transverse confinement is described by equally 
spaced harmonic-oscillator states, whereas the oscillator 
treatment is not applicable to the axial dimension due 
to the presence of the optical lattice. Our discussion 
hereafter is based on the tight-binding approximation 
that only the ground band is accessible to the system. 
This approximation is valid when the thermal energy 
of the atoms is much less than the first band gap of a 
deep lattice. For a ID uniform lattice, the eigen energy 
can be written as a function of quasimomentum q |15| . 
e (<z) = 5^ — 2 J cos(qd/h). Here, ui is the frequency of 
the local oscillation at each lattice well, while J is the 
tunneling energy due to the hopping to a nearest neigh- 
boring well, and it depends upon the lattice depth s in 
the following form 
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m 



J = —=E r s z/A exp (-2^5) 



(2) 



It should be noted that Eq. @ is valid only for deep 
lattices. At s = 11, for example, it overestimates J by 
approximately 18%. For the combined trap, it is a rea- 
sonable assumption that Eq. ([2]) remains valid as long 
as the trapping frequency uj z of the weak axial confine- 
ment is much smaller than the tunneling rate J/h. We 
are thus able to use a constant J over the entire lattice 
system at a given lattice depth. For simplicity the energy 
spectrum corresponding to the combined confinement in 
the axial direction is approximated by the semiclassical 
energy, e(p z ) + ^moj 2 z 2 , in the z-p z phase space, where 
p z is the quasimomentum in the ground band. Now we 
are able to write the total energy spectrum in an explicit 
form 



,{z,Pz) =f^±(n x + n y + 1) + ^mujlz 2 



1 



(3) 



huj z — 2J cos(p z d/h), 



where {n x ,n y } are non-negative integers. 

For a semiclassical description of this system, we treat 
the harmonic trap semiclassically while treating the op- 
tical lattice quantum mechanically. Such a picture cor- 
responds to a density distribution of the thermal cloud: 



E / ^:F(p z ,z)Md\%M\ 2 



(4) 



where 



V(x,y,z) 



^moj'iix 2 +y 2 ) + ^muj 2 z 2 + V sin 2 (kz). (1) 



exp[ 



M)] 
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and 



M/2 



>M 



w(z — Id) exp (ip z z/h) . 



l=-M/2 



Here, $ Pz is the normalized wave function of a uniform 
optical lattice system with an extension of M lattice sites, 
and w{z — Id) is the Wannier wave function. The total 
number of thermal atoms is then written as 



N, 



th 



E 



n(z)dz 

dp z dz 



2nh 



F(p z ,z)Md\% z {*t 



(5) 



In the tight-binding limit, w(z) is well localized within a 
single lattice site. In contrast, F(p z ,z) is a slowly vary- 
ing function of z. Therefore, $ Pi (z) in Eq. ([5]) can be 
integrated out. This results in a new integrand expressed 
as a summation of discrete F(p z , z — ld)d, which in turn 
can be approximated as an integral over z. By doing so, 
one gets 



N th = 



E 



dp z dz 
2ttH 



F( Pz ,z). 



(6) 



Below a critical temperature T c the chemical potential 
p, of the Bose gas reaches the bottom of the ground band 
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-huj z — 2 J, 



while the lowest state with p z — becomes macroscopi- 
cally populated which corresponds to the onset of Bose- 
Einstcin condensation. The condensed atoms exhibit 
macroscopic superfluidity, whereas all other atoms be- 
yond the lowest state form a so-called normal gas. Since 
the condensate is actually a quantum fluid, we use "su- 
perfluid" just as a synonym of BEC. The atomic number 
of the normal gas is given by the sum, 



1 



dp z dz 



exp[j3(s nxn - p c )} - 1 2irh 



N nc = N-N c = 

(7) 

where N is the total number of the atoms, N c the atomic 
number of the condensed part, j3 = X/ksT, and k B the 
Boltzmann constant. The integrand can be expanded in 
powers of the exponential term using the formula l/(e x — 
1) = Y^n=i e ~ nx ■ Moreover, the sum over n x , n y can be 
replaced by an integral if the atomic number TV is large. 
Performing the integration over n x , n y as well as that 
over the coordinate z, one gets 



N- N„ 



1 



oo ^ 
^E , 7 2 



(PhLOj_) 2 ^ n 2 

v / n—l 

dp. 



2n 



1/2 



(8) 



^ exp[-n/3(2 J - 2Jcos(p,d/»))], 



and the right side is a function of temperature T. Ap- 
parently, Eq. ([HJ is suitable for numerical calculation of 
the atomic number in the normal gas since the integra- 
tion can be simply replaced by a summation over the p z 
region. 

By imposing that N c = at the transition, Eq. (jHJ 
determines a critical temperature T c for a given N. Ap- 
parently, to obtain the value of T c , one needs to carry out 
numerical computations based on Eq. flSJ. Nevertheless, 
we can derive an analytic expression of T c in a limiting 
case. When the temperature of the Bose gas is so low 
that most atoms occupy the states in the vicinity of the 
bottom of the ground band, the relation p z <C h/d holds, 
and the cosine function in Eq. ((8]) can be expanded to the 
order of p 2 z . With the p z -dependent function integrated 
out, one has 

N = (k B T c /mf (m*/m) 1/2 C(3), 

where ZJ — (u^u^) 1 / 3 is the geometric average of the 
trapping frequencies, m* = h 2 /2Jd 2 the effective mass 
of the atom, and C( a ) = S^°=i 1/^" the Riemann zeta 
function. Finally, one gets 



k B T c = 0MhiJN 1/3 (m/m*) 1/6 ■ 



(9) 



which can be used as an analytic estimation of the critical 
temperature. 
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FIG. 1: Critical temperature T c versus the total number 
of S7 Rb atoms. The solid curve and the dashed curve are 
obtained from the numerical calculation of Eq. (JSj) and 
the analytical approximation of T c (Eq. ©), respectively. 
The dotted line gives the full numerical result by diagonal- 
izing the single-particle Hamiltonian. The lattice parameters 
are d = 400 nm and s = 11.2 E r . The trapping frequen- 
cies of the harmonic potential are uj± = 2ir x 83.7 Hz and 
u) z — 2ir x 7.63 Hz, respectively. 

We recall that an ideal Bose gas trapped in a 3D har- 
monic potential undergoes the phase transition to Bosc- 
Einstein condensation at a temperature [l6j /cbT c = 
0.94SoJiV 1 / 3 . Comparing this expression with Eq. ([9]), 
one can see that T c is changed by a factor of (m/m*) 1 ^ 6 
due to the presence of the ID lattice. Since m* is always 
larger than m Q under tight-binding approximation, the 
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combined trap T c is actually reduced compared to the 
case without lattice. Note also that a homogeneous 3D 
lattice system has a reduced T c as well [4|, but with a 
reducing factor ^Jm/m* instead. 

We have calculated the critical temperature T c for a 
87 Rb gas in the combined trap (see Fig. Q}. The trap 
parameters are intentionally chosen to match our ex- 
periment which will be described in the later section. 
The numerically calculated T c is displayed by the solid 
curve, while the dashed curve is the analytic T c calcu- 
lated according to Eq.©. The discrepancy between the 
two curves becomes larger as the atom number N is in- 
creased, showing that the accuracy of the analytic esti- 
mation becomes worse for larger N . We thus use only the 
numerically calculated T c in the following computations. 



tion. By diagonalizing the matrix (i H z j), one can get 

directly the energy spectrum e Uz . Furthermore, with the 
following formula, 



N : 



E 



i 



n x ,n v ,n z 



-n)/k B T 



1 



(12) 



we give the full numerical results of the critical tempera- 
ture and condensate fraction by the dotted lines in Figs. 
[1] and [H respectively. Clearly, our semiclassical treat- 
ment agrees with the full numerical method, and proves 
to be reliable. Furthermore, it offers a convenient way 
to analyze the spatial distribution of a confined atomic 
cloud, which in turn simplifies the calculation of interfer- 
ence patterns. 




FIG. 2: Condensate fraction as a function of T/T c . Solid 
line, numerical results based on Eq. ([8]). Dotted line is the 
full numerical result based on diagonalizing the single-particle 
Hamiltonian. The dash dot line is 1 - (T/T c ) 3 for T < T c . 
The parameters of the combined trap are the same as in Fig. 
[TJ and the atomic number is TV = 5 x 10 4 , corresponding to 
T c = 47.9 nK. 

It is well known that an ideal Bose gas in a 3D har- 
monic potential shows a T dependence of the condensate 
fraction as N c ~ 1 - (T/T c ) 3 for T < T c . We have also 
calculated the condensate fraction for our combined trap 
system with 5 x 10 4 atoms, as shown by the solid line in 
Fig. [U It displays a noticeable deviation from the curve 
of 1 — (T/T c ) 3 , but fits well to the characteristic shape, 
1 - (T/T C ) Q , with a = 2.679. 

To justify our analytical approximation, we also cal- 
culate the critical temperature and condensate fraction 
based on the diagonalization of the single-particle Hamil- 
tonian. The energy spectrum of the system is written as 

£ nx n v n z = Tw± (n x + n y + 1) + e„ a . (10) 

s nz can be obtained numerically from the following 
single-particle Hamiltonian along the z direction, 

H z = - - ^ (oJoj + «] ) + E £ '^ 1 ■ ( 1 1 ) 

(id) i 

Here £j describes an energy offset at each lattice site due 
to the presence of the harmonic trap along the z direc- 



III. INTERFERENCE PEAKS 

When the combined trap is suddenly shut off at the 
moment t = 0, the Bose gas starts to expand freely. Af- 
ter a time of flight r, the expanded wavepackets initially 
localized in single lattice wells overlap with each other, 
forming a 3D density distribution. In the following cal- 
culation, the x and y dependence of the atomic density 
will be integrated out so as to obtain a density profile 
along the z direction only. This is convenient for making 
a comparison with the experimental results. Usually, an 
absorption image is used to record the column density 
profile of a released atomic cloud. Supposing that the 
probe laser beam is applied along the direction of the x 
axis, the density profile along z can be easily obtained by 
integrating the column density over the y dimension. 

A. Normal gas 

In the combined trap, normal gas atoms are distributed 
over the transverse harmonic modes labeled by a positive 
quantum number q = n x + n y . For a given q, there are 
5 + 1 degenerate states, and we hereafter call them sub- 
states. The summation over n x and n y in the previously 
mentioned equations is thus equivalent to ^2 q {q + 1) ■ • • . 
From Eq. (J7]), one sees that the substates belonging to 
the same q number have identical populations and spa- 
tial distribution along the z direction. Due to optical 
lattice potential, atoms in a single substate are further 
distributed over the Bloch states with different quasimo- 
mentum p z with p z /h G (—ir/d,ir/d). Each p z compo- 
nent can be treated semiclassically where the influence of 
the optical lattice is given by a quantum wave packet de- 
scription, while the influence of the harmonic trap along 
the z direction is treated semiclassically. In such a pic- 
ture, the single-particle wave function of a p z component 
at t — takes the following form: 

*jL (* = 0) = a M z ~ ld ) exp(iPzz/K), (13) 
l 
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where (a*) 2 denotes the probability for a particle roughly 
located in the Ith lattice site for the transverse harmonic 
mode q. 

Equation ([7} shows that the atomic density of a sub- 
state with p z has an envelope as 



n(z) 



Ap z 



1 



2nK e i8[gfiwj_+2J(l-cos(p 8 d/ft))+imw2z 2 ] _ y 



(14) 



where Ap z denotes a small interval of p z 
number in the Ith lattice site is then 



ni 



dAp z 



2irh P l3[qtiu±+2J(l~cos(pzd/h)) + ±mu>ld 2 l 2 } 



The atom 



(15) 



Therefore, a\ is simply given by (a/) 2 = ni/N q , with 
N q = '^2 l ni being the total atom number of the p z com- 
ponent in the substate of interest. In principle, ol\ should 
be determined by solving the Schrodinger equation of H z . 
However, as shown lately, the thermal average of in-trap 
density written in terms of |af | 2 is matched to the ex- 
pression obtained by semiclassical approximation, hence 
within semiclassical approximation \af\ 2 can be identified 
to (c^) 2 =ni/N q . 

In the tight-binding limit w(z) can be well 
approximated by a Gaussian wave packet 
(ttu 2 )^ 1 ^ 4 exp(— z 2 /2a 2 ), where a = y/h/muj z is 
the oscillator length. After the free expansion over a 
time of r, the single-particle wave function of the atoms 
with p z is written as 



W (t = r) 



i 



K{z, z' , t)w(z' — Id) exp(ip z z' /h)dz' . 



(16) 



Here, K(z,z',t) = (z\ exp(— iHt jK)\z') is the propaga- 
tor, with H the Hamiltonian governing the expansion 
process. If the interatomic interaction is neglected, H 
contains only the kinetic energy, say, H — P 2 /2m. In 
this case, it is straightforward to get 



K(z,z',t) = 



',2nhT 



exp 



2Ht 



{z-z'f 



(17) 



For simplicity of expression and calculation, we will 
use dimensionless units for the length in z, the quasimo- 
mentum p z and the time t by the replacement z — > zd 
(and hence a — > ad), p z — > p z h/d, and t — > t(2md 2 /ti). 
Inserting Eq. (fTTj) into Eq. (fTS]). and working out the 
integration over z' , one gets 



*l(t = T)=Aj2a 9 l B(p z ,l,z), 



(18) 



with A and B given by 

A = n- 1 ^(a + i2T/a)- 1 ^ 2 , 



B( Pz ,l,z) 
exp 



-2tg 2 p 2 + 2a 2 zp z + i[4lp z T + (z - If 
2(tr 4 + 4r 2 )/(2r + ia 2 ) 



(a) I 




-1800 -1200 -600 600 1200 1800 



z/d 




FIG. 3: (a) The solid lines are the calculated atomic distribu- 
tion of the normal gas ( 87 Rb) after 30 ms of time of flight. The 
total atom number iV = 5.9 x 10 4 , and the trap parameters 
are the same as in Fig. [T] corresponding to T c — 51.1 nK. For 
the top two curves, the temperature T/T c is 0.98 and 0.95, 
respectively. The other nine curves, from top to bottom, are 
for T/T c ranging from 0.9 to 0.1 with a step of 0.1. Vertical 
scales of the curves have been adjusted so that the central 
peaks have roughly the same height, (b) Solid circles repre- 
sent the visibility of the side peaks in (a). The solid curve 
connecting the points is added to guide the eye. Open circles 
show the normal gas fraction for the given total atom number. 



Taking into consideration all transverse modes and all 
p z components, one gets the atom density after the time 
of flight, 



,( Z )^^( ? + l)JV,|fJJt = r) 



9=1 Pz 



9=1 Pz 



22^B(p z ,i,z) 



(19) 



In the numerical calculations of Eq. (|19[) . the summa- 
tion over transverse modes is cutoff at q = 200, while 
n is cutoff at 30 and lattice number I at ±350. These 
cutoff numbers are chosen to assure a high accuracy bet- 
ter than 0.2% in the calculations of atom numbers. By 
setting Ap z to 0.05-7T, the whole p z range is divided into 
40 intervals. This step size of p z has the order of h/dM, 
where M ~ 100 is the typical spatial extent in the z di- 
rection. We have also checked that the calculated results 
have almost no change when further reducing Ap z . The 
step size in z is set to be 18d (7.2 /im), comparable to the 
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pixel size of 9 /im in our experiment. 

In Fig. [3Ja) we show the numerical results of a nor- 
mal gas of 87 Rb atoms after 30 ms of time of flight. It 
is obvious that the normal gas leads to three peaks in 
the atomic distribution along the z dimension. However, 
these peaks are not sharp, and, close to the T c , the side 
peaks are not even well resolved. In contrast, a normal 
gas initially trapped in a 3D homogenous lattice system 
gives rise to much sharper peaks Q. 

We now define a visibility for the side peaks as in Ref. 
171: 



riA - n B 
n A - n B 



(20) 



where tia is the atomic density of the side peak, and ns is 
atomic density at the minimum between the central peak 
and the side peak. From Fig. [3]Jb) it can be seen that the 
visibility v is well below 1 at a considerable fraction of 
the normal gas. Although v can reach 0.8 at a very low 
temperature of T = 0.1T C ~ 5nK, only about 0.1% of 
the atoms remain in the normal gas state while all other 
atoms are condensed. Actually, it is hard to detect such a 
small number of atoms using the conventional absorption 
imaging method. 



B. Condensed gas 

Unlike the normal gas, Bose-condensed atoms in the 
combined trap occupy the lowest state of q = and pile to 
a small quasimomentum interval of p z = 0. Nevertheless, 
the normal gas propagator holds also for the condensed 
gas. In analogy with the calculations for the normal gas, 
one can derive the single-particle wave function of the 
condensed gas after the time of flight, 



T)=A(T)^a k B (k,z), (21) 



where k denotes the kth lattice well, a 2 is the probability 
of an atom staying in the fcth well, and 



Bo{k, z) = exp 



(i2t - a 2 ){z - kf 
2(a 4 + 4r 2 ) 



Note that, in the two formulas above, r, z and a are in 
their dimensionless form. 

In the tight-binding limit, condensed atoms in the com- 
bined trap form an array of subcondensates along the z 
axis. Each subcondensate is a 2D quantum gas in nature, 
and its density distribution in the radial dimensions is de- 
scribed by a Thomas- Fermi profile [18j . The local chem- 
ical potential associated with the repulsive interaction of 
the atoms depends upon the average atom number in the 
following form [19(: 



' gmu^Nk 
V2tt 3 /V ' 



(22) 



Here, Nk is the average atom number in the fcth lattice 
well, and g = Aith 2 a/m is the interaction parameter with 
a the s-wave scattering length. The sum of ^), oc and the 
external harmonic potential (l/2)muj 2 z 2 is just the chem- 
ical potential which should remain invariant throughout 
the entire condensed gas at equilibrium. Accordingly, Nk 
is given by 



N k = (15JV c /16*CAf)(l 



^ 2 Am) 



(23) 



where &m labels the outermost lattice w ell p opulated 
with condensed atoms, and it is written as [181 ] 



h 2 



2m fUN c a d 
nu> 2 d 2 \ 8i/n a/ lo a 



2/5 



(24) 



If we neglect the mean- field interaction of the condensed 
gas during the free expansion time, the density distri- 
bution after the TOF can be directly derived from the 
coherent superposition of the expanded subcondensates: 



n c (z) = |A(r)p 



\/Nk~B (k,z) 



k=— fe_j 



(25) 



For the condensed gas before expansion, due to the 
high atomic density, /i)? c /^ is in the order of several hun- 
dred Hertz. Although the atomic density drops quickly 
during the TOF, we still expect that mean-field interac- 
tion might lead to a considerable change in the coherence 
property of the expanding atomic clouds. For simplicity, 
we only consider the mean-field interaction within the 
single expanding subcondensates, and neglect the inter- 
action between them. At the beginning of the TOF, the 
combined trap is suddenly turned off. Therefore, the to- 
tal energy of the fcth subcondensate includes only the 
mean-field energy at this moment, that is, 

E k = E int = (l/2)U k N 2 . 

Here, Uk = g J |$fc(r, z)| 4 dzdr is the on-site interaction 
matrix element of the fcth subcondensate when confined 
in the combined trap. Using the analytic form of <l>fc(r, z), 
we can get the expression of Uk in terms of the trap 
parameters: 



(26) 



Its dependence on Nk is due to the fact that the atomic 
number affects the Thomas-Fermi radius of the radial 
wave-function. At later times, the total energy Ek 
remains constant despite the fact that the interaction 
energy is being converted into kinetic energy. Then, 
the corresponding chemical potential is simply given by 
fik = dEk/dNk = a4°7 2 - Over the total time of flight, 
the fcth subcondensate acquires an additional phase fac- 
tor exp(— inkT/h). Consequently, we can write the den- 
sity distribution at the end of TOF by just inserting this 
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phase factor into Eq. ([25|) : 



n c {z) = \A{r)f 



£ 

k— — kj\, 



'Nkexp(-ifj,kT/h)B Q (k,z) 



(27) 

From this equation one sees that fj>k will affect n c (z) by 
its nonuniformity. In Fig. EJa), we plot two curves cal- 
culated, respectively, with Eqs. (|25j) and (f27| for a con- 
densed gas of 2.5 x I0 4 atoms. When the mean- field in- 
teraction is included, all three interference peaks are sig- 
nificantly broadened by about a factor of two. Since the 
mean-field interaction is nonnegligible, all the theoretical 
interference patterns mentioned hereafter are computed 
by Eq. (@7|). 



1800 




-1800 -1200 -600 600 1200 1800 

z/d 



FIG. 4: Atomic distribution of the condensed gases after 
30 ms of time of flight. The total atom number and the trap 
parameters are the same as in Fig. [3] (a) The dashed (solid) 
line is the result for T — 0.7T C with (without) a consideration 
of the mean-field interaction during TOF. (b) These curves 
are calculated from Eq. (|27|) . From the top curve to the bot- 
tom one, the temperature is decreased in sequence, with the 
same values as in Fig. [3] Vertical scales of the curves have 
been adjusted so that the central peaks have the same height. 

In Fig. 0Jb), we plot the calculated density distribu- 
tions of the condensed gases at different temperatures for 
a fixed total atomic number and a fixed lattice depth. 
With decreased temperature, the condensate contains 
more atoms, leading to wider interference peaks due to 
mean-field interactions. Unlike a normal gas, a con- 
densed gas always shows fully resolved interference peaks, 
with a high visibility very close to 100%. This charac- 
teristic behavior can be easily understood as the global 



coherence of condensed atoms in a superfiuid state. Addi- 
tionally, these peaks are considerably narrower than that 
of the corresponding normal gases, except the extreme 
cases of very low temperatures that the normal-gas atom 
number is very small and hardly detectable. When one 
measures the interference pattern of a mixture of the con- 
densed gas and a normal gas, one would see three narrow 
peaks riding on three broad peaks. This is the so-called 
"peak on a peak" structure which was first predicted for 
a homogeneous system [4, 5]. For an inhomogeneous sys- 
tem in the combined tap, the onset of the condensed gas 
is also characterized by the "peak on a peak" structure. 

On the other hand, if the condensed gas undergoes 
only a ballistic expansion during the TOF (no mean-field 
interaction) , the relative intensity of the side peaks with 
respect to the central one should obey a simple law [l8[ : 
Pi = exp(-4/v/i). We check the data in Fig. gjb) (Pi = 
0.303 for s — 11.2), and find that the side peak intensities 
agree well with Pi (to within 2 percent). It seems that 
the analytic expression of Pi is also valid in the case of the 
existence of mean-field interaction during the expansion 
time. 



IV. EXPERIMENT 

In experiment, we create a cigar-shaped 87 Rb conden- 
sate in the hyperfine state |P = 2, Mp — 2), confined in 
a conventional Quadrupolc Ioffc Configuration (QUIC) 
trap with an axial trapping frequency of uj z = 2tt x 18.7 
Hz and radial trapping frequency of lu±_ = 2it x 205 Hz. A 
nearly pure condensate contains approximately 2 x 10 5 
atoms. If the frequency of the rf knife for evaporation 
cooling is ramped down further, we can obtain a conden- 
sate with a lower temperature at the cost of decreased 
atomic numbers. Certainly, the temperature is hard to 
measure because there are almost no thermal atoms re- 
mained. Nevertheless, we are able to coarsely adjust 
the temperature of the cold atomic sample using the 
rf knife. After the evaporation cooling, the QUIC trap 
is adiabatically relaxed until the axial and radial trap- 
ping frequencies are lowered to ujj_ = 2~k x 83.7 Hz and 
lu z — 2tt x 7.63 Hz, respectively. Accordingly, the spatial 
extension of the condensate wave packet is increased by a 
factor of 2.45, so as to cover more lattice periods at later 
times. The optical lattice is formed by one retroreflected 
laser beam which is derived from a Tksapphire laser at a 
wavelength of A = 800 nm and focused to a 1/e 2 radius 
of 300 /im. It is applied to the condensate along the long 
axis, and it is ramped up to a given intensity over a time 
of 50 ms and held at this value for 10 ms. The sum of the 
optical lattice and the QUIC trap potential gives a com- 
bined trap in accord with Eq. ([TJ. The potential depth 
of the optical lattice is calibrated using the method of 
Kapitza-Dirac scattering [l7|. We then suddenly switch 
off the combined trap and allow the cold atomic sample 
to expand freely for a time of 30 ms. Finally, we take 
an absorption image of the expanded atomic gas using a 
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CCD camera, from which we can deduce both the total 
atom number and the atomic density distribution. 
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FIG. 5: Density distribution of two typical atomic samples 
of 87 Rb. Black curves are the measured linear density of the 
released 87 Rb atoms after a 30 ms TOF. Gray curves are the 
calculated results based on Eqs. (I19[) and (|27l) . where tem- 
perature is used as a fitting parameter. The dotted curves 
in (b) and (c) are the calculated results based on measured 
temperature values. Their vertical scales have been adjusted 
for comparison with the black curves, (a) The measured total 
atom number and lattice depth are 7V~5.3xl0 4 ands~5.6, 
respectively, corresponding to T c = 55.9 nK. T = 44.7 nK is 
assumed in the computation, (b) N ~ 5.9 x 10 4 , s ~ 11.2 and 
T c = 51.1 nK; T = 33 nK for the gray line and T = 37.4 nK for 
the dotted line, (c) N ~ 1.1 x 10 5 , s ~ 16.7, and T c = 63.0 nK; 
T = 55 nK for the gray line and T = 49.3 nK for the dotted 
line, (d) TV ~ 2 x 10 5 , s ~ 20, and T c = 76.5 nK; T = 73.4 nK 
for the gray line. 

The "peak on a peak" features of interference patterns 
were confirmed by the measured linear densities of ex- 
panded atomic clouds in many runs of experiments. The 
black curves in Figs. Efa)-[5]Jc) display three typical den- 
sity distributions which were obtained by integrating the 
pixels in each column of the corresponding absorption 
image. The calculated critical temperature T c is usu- 
ally in the order of several tens of nK. In contrast, dur- 
ing the evaporative cooling stage, the critical tempera- 
ture for condensation in the QUIC trap is much higher 
(~ 400 nK). Actually, the experiment reported in [201 ] 
showed clearly a significantly decreased critical temper- 
ature for a combined trap when compared with a purely 
magnetic trap. 

To test our theory, the temperature of the atomic sam- 
ple before expansion must be known. It can be deduced 
from the transverse distribution of the normal gas after 
the TOF, which should take a Gaussian profile due to 
the initially thermal occupation of the transverse modes. 
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FIG. 6: Density distributions in the transverse dimension for 
the same atomic cloud as in Fig. [5jb). The black curve is the 
integration along the axial direction of the central peak. The 
gray curve is a bimodal fitting of the black curve, which rep- 
resents a superposition of a Gaussian profile and a Thomas- 
Fermi profile (inverted parabola). The rms size of the Gaus- 
sian part gives a temperature of 37.4 nK. (Inset) The absorp- 
tion image of this atomic cloud. 



Since the condensate part should take a Thomas-Fermi 
profile in the transverse direction, a bimodal transverse 
distribution is expected for a released gas, just as shown 
in Fig. |Hl However, a measured temperature based on 
this method usually has a large uncertainty due to the 
following reasons. First, in the transverse direction, the 
spatial extent of the condensate is not very distinct from 
that of the normal gas, especially when the condensate 
fraction is large. Second, the normal gas density profile 
deviates from an ideal Gaussian shape, and exhibits a 
slight asymmetry that may arise from the misalignment 
between the lattice light and the magnetic trap. Third, 
the optical noise in the absorption ima ges also lowers 
the fitting accuracy. As pointed out in [2l|, an atomic 
sample can be significantly heated or cooled when adi- 
abatically loaded to an optical lattice. Yet, to date, we 
have no alternative methods for accurate measurement 
of the temperature of an atomic sample confined in a 
lattice system. We have to treat temperature as a fit- 
ting parameter in the calculation, so that the calculated 
density distribution most closely reproduces the experi- 
mental curve. 

The gray curves in Fig. [5] are the calculated density 
distributions of cold rubidium gases. Figure [5£ a) is a 
case with a smaller atom number, and the fitting curve 
agrees fairly well with the experimental data. Figure 
[5Jb) displays the interference pattern of another atomic 
sample initially confined in a deeper lattice, and the side 
peaks are more prominent. As a comparison, the atomic 
sample in Fig. [5jc) contains more atoms and the lattice 
depth was further increased. Accordingly, the calculated 
T c shifts up to 63 nK. Since the normal gas density is 
increased, the feature of "peak on a peak" is more pro- 
nounced. In (b) and (c), we also plot the density distribu- 
tions calculated using the measured temperature values. 
The larger deviation from the measured density profiles 
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should be attributed to the inaccuracy of temperature. 

Despite the overall agreement between the theoretical 
and experimental curves in Fig. [5j there are still no- 
ticeable discrepancies. As the lattice depth increases, 
theoretical normal gas peaks become broader than the 
measured density profiles. As shown in Fig. EJd), the 
gray line does not match the black curve, particularly 
at the wings of the normal gas. For this atomic sam- 
ple, the temperature is close to T c , and the condensate 
peaks are hence very small. The mismatch between the 
theory and experiment indicates that our model is not 
valid for very deep optical lattices. This can be easily 
understood by comparing the tunneling energy J to the 
energy offset between adjacent lattice sites induced by 
harmonic-potential. For a cloud with an extension of 
Im<I, this energy offset is muj^d 2 lM at the site of Im- A 
typical value of Im = 200 corresponds to an energy offset 
of 2nh x 16 Hz, whereas J decreases with increased lat- 
tice depth. For s = 20 as in Fig. EJd), J ~ 2nh x 10 Hz. 
When J gets smaller than the energy offset between adja- 
cent lattice sites, normal gas atoms are essentially local- 
ized and the semiclassical analysis breaks down. We have 
also calculated the T c 's for situations of s > 20, using the 
diagonalization method and semiclassical approximation, 
respectively. We do find significant discrepancy between 
the predictions of these two methods. In such situations, 
an atomic cloud should be treated as a chain of separate 
condensates, where the loss of condensate is interpreted 
as the loss of well-to- well phase coherence [13] ■ At a 
depth level of s ~ 30, we observed a completely disap- 
pearance of interference peaks. 

On the other hand, in our model, the interatomic in- 
teractions during TOF are taken into account for the 
subcondensates individually. This amounts to neglecting 
the variation of the wavefunction modulus induced by the 
interactions between subcondensates. The fluctuations 
of the wavefunction modulus are related to the relative 
phase of the subcondensates, and hence affect their phase 
coherence, leading to variations of the interference peak 
of the condensates. Evidently, the calculated condensate 
peaks are slightly wider than the measured ones (see Fig. 
[5]). Perhaps, the neglected interactions are favorable for 
establishing a uniform phase which partially cancels the 
phase nonuniformity discussed in Sec. 1111 Bl Since we 
have not found a simple model to account for it, this 
effect will not be discussed in detail in this paper. 



V. CONCLUSION 

We have performed a study, both theoretically and ex- 
perimentally, on the phase transition to macroscopic su- 
perfluidity for a Bose gas confined in a combined trap 
formed by a harmonic potential and an optical lattice. 
We have mainly investigated the interference patterns 
of the Bose gases below the critical temperature. By 
using a semiclassical energy spectrum and tight-binding 
approximation, we have derived an analytical approxi- 
mation of the critical temperature which is applicable 
to an atomic gas residing in the vicinity of the bottom 
of the ground Bloch band. Furthermore, the interfer- 
ence patterns of the normal gas and the condensed gas 
were analyzed separately. We have derived the analyti- 
cal expressions of the atomic density distribution for the 
released normal gas and condensed gas which has ex- 
perienced a free expansion over a time of flight. Our 
calculation results show that a condensed gas is charac- 
terized by fully-resolved narrow interference peaks while 
a normal gas forms broad interference peaks with lower 
visibility. For comparison, we have performed a prelimi- 
nary experiment using Bose-Einstein condensates of 87 Rb 
atoms. The combined trap system was implemented by 
applying a ID optical lattice to a magnetically trapped 
condensate. The measured interference pattern agrees es- 
sentially with our theoretical prediction, exhibiting "peak 
on a peak" structures associated with the onset of con- 
densed gases. Thus, both the theoretical and experimen- 
tal results confirm that the "peak on a peak" structure 
is indeed a signature of macroscopic superfluidity in our 
inhomogeneous lattice system. 
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